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Abstract 

We integrate numerically the Kardar-Parisi-Zhang (KPZ) equation in 1 + 1 and 2 + 1 dimensions 
using an Euler discretization scheme and the replacement of (V/i) 2 by exponentially decreasing 
functions of that quantity to suppress instabilities. When applied to the equation in 1 + 1 dimen- 
sions, the method of instability control provides values of scaling amplitudes consistent with exactly 
known results, in contrast to the deviations generated by the original scheme. In 2 + 1 dimensions, 
we spanned a range of the model parameters where transients with Edwards- Wilkinson or random 
growth are not observed, in box sizes 8 < L < 128. We obtain roughness exponent 0.37 < a < 0.40 
and steady state height distributions with skewness S = 0.25 ± 0.01 and kurtosis Q = 0.15 db 0.1. 
These estimates are obtained after extrapolations to the large L limit, which is necessary due to 
significant finite-size effects in the estimates of effective exponents and height distributions. On the 
other hand, the steady state roughness distributions show weak scaling corrections and evidence of 
stretched exponentials tails. These results confirm previous estimates from lattice models, showing 
their reliability as representatives of the KPZ class. 
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I. INTRODUCTION 



Nearly two decades ago, the Kardar-Parisi-Zhang (KPZ) equation 

— = U2 ^ h + X 2 (Vh) 2 + r ] (x,t). (1) 

was proposed as a hydrodynamic description of interface growth In Eq. ^TJ), /i is the 
interface height at position x and time t, the linear term represents the effect of surface 
tension, the nonlinear term accounts for an excess velocity due to local slopes and r] is a 
Gaussian noise with zero mean and co- variance (r/ (x, t) r] \x' ', t'^j) = D5 d (x — x'^j 5 (t — £'), 
where D is constant and d is the dimension of the substrate. Several applications of the 
KPZ equation are catalogued in Refs. and recent examples in d = 1 and d = 2 are 

presented in Refs. {5, 6]. These applications and the intrinsic interest as a non-equilibrium 
statistical mechanical model motivated a intense theoretical study of its properties and of 
properties of lattice models in the KPZ class, i. e. models which obey the KPZ equation in 
the continuum limit (long times, large sizes). 

Many properties of KPZ systems in d — 1 are known because the steady state height 
distribution for periodic boundaries is the same of the Edwards- Wilkinson equation (the 

n 

case A = - EW) [7[. This result and the Galilean invariance property provide the exact 
values of the scaling exponents of the average surface roughness 111 M. The full height 

Tin 

distributions for other boundary conditions in d =^_1 are also known [5J, [8[ , 



but controversies 



on the universality of correlation functions exist 



9| and are subject of current work 10] 



On the other hand, a small number of exact results are known in the most important 



case for applications to real systems, which is d = 2 
in d > 2 from some analytical approaches 
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Exponents estimates were obtained 



deviate from accurate numerical results of lattice models 
distributions were also calculated numerically 



15 
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but their predictions usually 



18 



16l . Il7j . Height and roughness 
T9I . [20! because they may 



be useful for comparison with real systems data and for additional tests of the analytical 
predictions. However, the accurate numerical data currently available were obtained only 
from two or three lattice models, such as the restricted solid-on-solid (RSOS) model [21] . 
because those works aim at reducing scaling corrections to improve the accuracy of the final 
estimates. Consequently, no systematic variation of the parameters of the KPZ equation 
can be performed in such works (although the parameters of the KPZ equation associated 



to each lattice model may be determined by inverse methods 



22 



23|]). 



It is certainly desirable that the universality of the above mentioned quantities is also 



tested with the KPZ equation itself with a suitable variation of its coefficients, 
integration of the KPZ equation was already performed by several authors 



29. 
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ndeed. 



25 



the 



26 



30 



n 



3l| . However, they usually focus on the calculation of scaling exponents (which 



typically have lower accuracy than the discrete models data). Giada et al [29J discussed the 
relevance of other quantities to characterize KPZ scaling, such as the skewness of height 
distributions, but they did not determine their universal values. The most recent work on 
the subject seems to be that of Ma et al {30]], which suggests exponent values very different 
from the previous ones. 

The aim of this work is to fill that gap by analyzing numerical results for height and 
roughness distributions in the steady state of the KPZ equation with several sets of coeffi- 
cients. Estimates of roughness exponents will also be provided here. Our results confirm the 
universality of these quantities, with the values previously suggested by lattice model simu- 
lation. Of particular relevance is to confirm that the roughness distribution has a stretched 
exponential tail, which reflects the non-Gaussian behavior of the interface. We also show that 
the calculation of reliable roughness exponents and of dimensionless amplitude ratios char- 
acterizing the height distribution has to account for the presence of finite-size corrections, 
which are much smaller in the scaling of roughness distributions. These features resemble 
those found in lattice models and show that the finite-size corrections are not artifacts of 
those models. 

We will adopt an Euler scheme for integration. The main problem of using a simple 
version of this method 27] is the onset of instabilities in the growing interface at long times, 



32 



331 ] . However, we 



CD) by an exponentially 



32 



331 ]. Another previ- 



as indicated by the divergence of the height at a certain position 
will control this instability by replacing the square gradient in Eq. 
decreasing function of this quantity, as suggested by Dasgupta et al 
ously reported problem of the simple Euler scheme is the anomalous value of the amplitude 
of steady state roughness scaling in d — 1 34|. However, we will show that this anomaly 
also disappears with the introduction of the exponentially decreasing nonlinear term, with- 
out needing special discretization schemes. The use of a framework which avoids different 
types of anomaly in the integration of the KPZ equation is certainly relevant and is a subject 



of current interest even for studies in d 



3l|. 



The rest of this work is organized as follows. In Sec. II we present the quantities of 



3 



interest of this work. In Sec. Ill we present the integration scheme and show results for the 
one-dimensional case, where exact solutions are known. In Sec. IV we present the results in 
d = 2. In Sec. V we summarize our results and present our conclusions. 



II. THE QUANTITIES OF INTEREST 

The simplest quantitative characteristic of a given interface is its roughness [2|, also called 
the interface width, defined as the rms fluctuation of the height around the average position. 
The squared roughness, w 2 = h 2 — h 2 , is usually averaged over different configurations, and 
its scaling on time and length is used to describe non-equilibrium growth processes. For 
short times, the average roughness scales as 

M ~ Bt 2 ?, (2) 

where f3 called the growth exponent and B is a scaling amplitude. For long times, in a finite 
system of size L, a steady state is attained, with the average width saturating at 

Msat « AL 2 «, (3) 

where a is called the roughness exponent and A is another scaling amplitude. 

The exponents a and (3, as well as the dynamical exponent z — a//3, are the quantities 
most frequently used to characterize a given universality class of growth. For the KPZ class, 
Galilean invariance leads to the additional exact relation a + z = 2 

A better description of the interface is provided by the full height distribution, measured 
relatively to the average height. The moments of the steady state height distribution, 



W n = ((h-h) ), (4) 



may be used to characterize it. Numerical works usually consider some dimensionless am- 
plitude ratios for this purpose Q, [l^, particularly if finite-size effects in the sea 
distributions are found and extrapolations to infinite system size are necessary 
lowest order ratios are the skewness 



ed height 
Q. The 



which is related to the asymmetry of the distribution, and the kurtosis 



(5) 



Q = ™\ - 3, (6) 



which is related to the weight of the tails of the distribution relatively to a Gaussian. 

Recent works suggest that the statistics of global quantities may be more useful for 
characterizing an interface growth problem. The main quantity of interest is certainly the 



squared roughness 



35 



36 



371 ]. whose probability of being in the range [1V2, w 2 + dw 2 ] will 



be denoted by P w (w 2 ) dw 2 . The probability density P w is expected to scale as 

Pw(w , ) = UhlzM\ t (7) 



where a = y (W2 2 ) — (w 2 ) 2 is the root-mean-square deviation. Compared to other scaling 
forms for P w (w 2 ), Eq. ([7]) has the advantage of being less sensitive to finite-size effects and, 
consequently, more useful in data collapse work [38]. Anyway, comparison of dimension- 
less ratios such as the skewness and the kurtosis of roughness distributions are important 
quantitative tests. 



III. INTEGRATION METHOD AND RESULTS IN d = 1 



The usual discretization of Eq. (TTJ) follows the lines of Ref. [27J , which gives in d — 1 



h(t + At) - h(t) 



At 



1 



(Ax) 



■{v [h (x - 1) - 2h (x) + h(x + l)} + -X[h (x + l)-h(x- l)] z } 



W12A1R (t) 



(8) 

In Eq. flS]), a = ^2D/{Ax) d and R is a random number taken from an uniform distribution 
in the interval [—0.5, +0.5]. In 2 + 1 dimensions, similar contributions of the y direction to 
the laplacian and to the square gradient are added to the right side of Eq. ([8]). 



The spatial step Ax = 1 can be used without loss of genera 



ity 



since decreasing Ax 



would be equivalent to decreasing the parameter A in Eq. (TTJ) [26|, |27j, |34[ . Here, the lattice 
size L used in the numerical integration has maximum values L = 128 both in d — 1 and 
d = 2. On the other hand, At has to be sufficiently small to provide accurate results. One 
ensures that a certain value of At is suitable for a certain set of parameters by verifying 

ts. As usual, we adopt v, a and 



27]. 



that further decreasing of its value does not change the resu. 
g = X 2 D / v 3 as free parameters in the discretized equation 

Although Eq. ([S]) provides reliable estimates of growth and roughness exponents in d = 1 
with relatively small system sizes, this discretization procedure has some problems. First, 



Lam and Shin 



34( showed that it provides incorrect values of the scaling amplitude A in 



Eq. ([3]), which depends on v and D. The deviation from the expected exact value of A 
is observed even after extrapolation to L — > oo. The same authors proposed an improved 
discretization in 1 + 1 dimensions 39|, but it was based on particular properties of the 



KPZ equation in that dimension and, consequently, cannot be extended to the main case of 



interest here (2 + 1 dimensions). Secondly, Dasgupta et al [32j, |33| showed that the above 
discretization generates instabilities in the interface, in which isolated pillars or grooves grow 
in time on an otherwise flat interface. In d — 1, these instabilities typically appear for large 
values of g and for sufficiently long integration times and/or lattice sizes. Moreover, these 
instabilities are not consequence of unsuitably large time steps, but an intrinsic feature of 



the discretization of the KPZ and other nonlinear equations 331 ] . 



In order to solve the second problem (the instabilities), here we adopt the scheme proposed 



in Refs. 



32 



in which (V/i) 2 in the KPZ equation is replaced by / ((V/i) 2 ), where 

/(a.) = (l_e— )/c (9) 

with c being an adjustable parameter. This method avoids that large local height differ- 
ences lead to very large growth rates, which is the origin of the instabilities 33]. The new 
discretized equation is obtained along the same lines of Eq. (JBJ), i. e. the square gradient 
is estimated from the nearest neighbor height differences in all spatial directions and the 
corresponding value of / ((V/i) 2 ) is calculated from it. 

Notice that the replacement of (V/i) 2 by / ((V/i) 2 ) corresponds to the introduction of 
an infinite series of higher-order nonlinear terms in the KPZ equation. Their introduction 
does not change the scaling exponents and other universal quantities Q]. 

Figs, la and lb illustrate the growth of an instability when Eq. (jSJ) is used with v — 1, 
D — 1 and g = 48, in a one-dimensional interface with L = 500, when time increases from 
55.5 to 55.8 (At = 0.05 there). This instability disappears after the replacement of (Vh) 2 
by / ((^^) 2 ) with c = 1 and integration with the same parameters. It is important to 
stress that c cannot be very large, otherwise nonlinear effects become very weak and a long 
transient with EW scaling is found (similarly to what happens in lattice models - see e. g. 
Ref. 4o})- Anyway, in all our simulations in d = 1 and d = 2 using c not too small, no 
instability was observed in the growth regimes nor in the (very long) steady states. 

During the integration of Eq. (jSJ) with and without the instability control, we also 
analyzed the first problem mentioned above. Restricting the comparison to situations where 
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no instability is observed (small A and small lattice sizes), we measured the amplitude 



A' ee 12* (w 2 ) sat /L (10) 



(Eg 



34|- 



in the steady states. It is expected that A' — > 1 as L — > oo (a = 1/2 in d = 1) [4l| |- 

The finite-size estimates of A' are shown in Fig. 2. With the simple Euler method 
[8]), the numerical value of A' converges to a value close to 0.85, as reported in Ref. 
However, this discrepancy is also eliminated with the new discretization, i. e. with / (^(Vh) 2 ^ 
replacing (V/i) 2 . Fig. 2 clearly shows that the finite-size estimates of A' are consistent with 
an asymptotic value A' = 1 within small error bars. 

Thus, the method to control instabilities also solves another problem related to the dis- 
cretization of the KPZ equation, which is the incorrect estimation of scaling amplitudes in 
the regime where the original discretization [Eq. (jSJ)] seems to be stable. This advances 
over most previous works on the subject because, as far as we know, they analyzed those 
anomalies separately. The only exception seems to be a recent work which compared the 



original Euler scheme and pseudospectral methods [3j|, which shows that the latter avoids 
instabilities in d = 1 and provides the correct value of the amplitude A' [311] . One difference 
from the present approach is that here the KPZ equation was modified in real space. An- 
other important difference is that instability suppression with the pseudospectral method 
requires small time steps such as At ~ 10 -3 , while with At ~ 10~ 2 they have a non-negligible 

n 

chance to appear in d — 1 [3JJ]. On the other hand, here we obtained satisfactory results 
with At > 10~ 2 in d = 1 and d = 2. 



IV. RESULTS IN d = 2 

Here we consider four sets of values of the model parameters, with v = 0.5 and a = 
0.1 kept fixed and varying g = X 2 D/i> 3 , which corresponds to different intensities of the 
nonlinearity. Suitable values of the constant c were chosen to avoid instabilities, typically 
increasing with g. The values of the parameters used in each data set are shown in Table I, 
in increasing order of nonlinearity from A to D. 

Our aim is to span a reasonable region of the parameter space, but both very small 
and very large g are difficult to work with. Small nonlinearities must be avoided because 



the results would show long transients with EW scaling 



401 ] or with random growth. In 
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such cases, KPZ scaling would only be observed in very large box sizes, where it is difficult 
to generate a large number of independent steady state configurations due to the large 
saturation times (z « 1.6 [17]). On the other hand, working with very high g is not good 
because the large values of c introduce large irrelevant terms in the KPZ equation which 
play a role in the finite-size scaling, possibly showing crossovers from other dynamics. 

We worked with five different box sizes, ranging from L = 8 to L = 128. The time step 
here was At = 0.04, and Ax = Ay = 1, for all sets of the parameters. Maximum simulation 
times ranged from 10 3 for the smaller sizes to 10 for the largest one. These times are much 
larger than the saturation times, so that long steady state regimes were obtained in each 
realization. The total number of realizations for each set of model parameters was 2 x 10 3 
for 8 < L < 64 and 10 3 for L = 128. 

Finite-size estimates of the roughness exponents were calculated as 

m -l hN (L)/(w 2 ) (L/2)} 
a(L) = 2 ffi2 • (U) 

These effective exponents are expected to converge to the dominant, asymptotic exponent 
a as L — > oo. However, work on discrete models and with growth equations show that Eq. 
([3]) may have correction terms, which leads to an L-dependence of a(L). 

In Fig. 3 we show ce(L) versus 1/L for all sets of parameters of the KPZ equation. Since 
we were not able to obtain results in large box sizes in a reasonable computation time, the 
estimates of a are not so accurate as those calculated from lattice models Q]- F° r small 
nonlinearity (set A), we note a significant size dependence of thus the extrapolation 

for L —>■ oo has lower accuracy than that for sets B and C. The trend of the data for set D 
(high nonlinearity) is different from the other sets, probably due to a much more complex 
crossover to the asymptotic scaling. 

For the above reasons, our extrapolation of the effective exponents is mainly based on 
the behavior of the data sets B and C, which give 0.37 < a < 0.40 (we adopted the same 

sets) . These values agree with the best current 
17] . which suggests that those models actually 



3, r 



extrapolation methods of Ref. [17] for al 
estimate a ~ 0.39 from the lattice models 
work as representatives of the KPZ class. The trend o 
the simple rational guess a = 2/5 (proposed in Ref. 
discarded from the final error bar. 

The former estimates of a from integration of the KPZ equation were as low as 0.18 



1 data sets in Fig. 3 suggest that 
) is not valid, although it is not 



24] 



8 



and 0.24 



251 ]. but subsequent works provided estimates closer to 0.4 [26j, |28j. The exponent 



(3 = 0.24 obtained by Moser et al 



271 ] is consistent with the latter values (using a + z = 2) . 



A recent work with a pseudospectral method provided estimates ranging from 0.38 to 0.40 
(central estimates) from the scaling of different quantities. Thus, it is surprising that the 
most recent numerical solution of the KPZ equation |30j suggests that z = 2 and that the 
exact result a + z = 2 is not obeyed. We believe that this discrepancy is caused by the use 



of very low nonlinearities in Ref. 



301 ]. which leads to a long EW regime. In this case, KPZ 



scaling can only be detected in very large boxes and very long times. 

Now we turn to the analysis of the height distributions. Due to finite-size effects, we 
focus on the scaling of dimensionless amplitude ratios that characterize those distributions. 
Thus, in Figs. 4a and 4b we show the skewness and the kurtosis of the height distribution, 
respectively, as a function of 1/L. 

Estimates of S are accurate and consistent with an asymmetry in the distribution, with 
positive S meaning sharp peaks and flat valleys (for positive A). Since S = for EW growth, 
a small value of S is a signature of a crossover from EW to KPZ. This is consistent with the 
results for L = 8 and L = 16 in Fig. 4a, which show increasing 5* for increasing nonlinearity, 
with 5* ~ 0.05 for the smallest g (set A). 

Universality of S in the continuum limit is suggested by extrapolation of the results 
of all sets to 1/L — > 0. The intersections of at least two extrapolated values, with the 
corresponding error bars, lead to S = 0.25 ±0.01. The central estimate is slightly below the 
value \S\ = 0.26 ± 0.01 obtained in the discrete models [Ha, ll7(], suggesting that \S\ = 0.25 
up to two decimal places. Here we recall that negative A would lead to negative S, but with 
an universal \S\, as discussed in Ref. [17]. 

The estimates of the kurtosis Q are less accurate, mainly for the largest sizes, but extrap- 
olations indicate a universal positive Q. Intersections of at least two extrapolated values 
lead to Q = 0.15 ±0.1, which is also consistent with the estimate Q = 0.134 ±0.015 obtained 

HQ. 



in discrete models 



Now we discuss the roughness distribution scaling. In Fig. 5 we show the scaled distri- 
butions [using Eq. ([7])] for two sets of parameters (C and D) in box size L = 128 and the 



distribution for the RSOS model in lattice size L = 256 [20j. The excellent agreement of 
these curves in three decades of the scaled probability density illustrates the universality 
previously suggested by simulation of discrete models 20J . Comparison of results in L = 64 
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and L = 32 show that finite-size effects are very small. Quantitative evidence of the agree- 
ment is provided by the estimates of their skewness and kurtosis: averaging results for all 
sets, we obtain S = 1.73 ± 0.04 and Q = 5.6 ± 1.1, which must be compared with lattice 
models data S = 1.70 ± 0.02 and Q = 5.4 ± 0.3 2(|. Here, the largest deviations from the 
central values are provided by the sets A and D, similarly to the other quantities. However, 
those deviations are so slight that they cannot be detected by visual inspection of the dis- 
tributions (see e. g. the data for set D in Fig. 5). Overall averages of S and Q are fully 
consistent with universality of roughness distributions, and the above discussion reinforces 
the conclusion that they exibit much smaller corrections to scaling than other quantities. 

An important feature of the KPZ roughness distribution is the apparently stretched 
exponential tail, which is suggested in Fig. 5 by a small upward curvature in the right 
tail. In order to analyze the tails of our curves, we assume that \I/ = oPi (1V2) decays as 
\l/ (x) ~ exp (—Ax 1 ), where x = W2 ~j W2 ^ [ see Eq. (JTJ)]. Thus, estimates of the exponent 7 
can be obtained from 

, [ In (*(*)) 

TW = n r ../ >■■ XTT> ( 12 ) 



ln(*(s-A)) 



In [x/ (x - A)] ' 
with constant A. 

In Fig. 6 we show 7 (x) versus 1/x 2 for three sets of parameters (B, C and D) and L = 64, 
using A = 4 in Eq. (1121) . This box size was used because fluctuations in the tails of the 
distributions for L = 128 are much larger. The trend of the data as x — > 00 suggests that 
the tail of the roughness distribution is an stretched exponential with an exponent between 
7 = 0.7 and 7 = 0.9. This is consistent with results of lattice models, which give 7 rs 0.8. 



V. CONCLUSION 



We solved numerically the KPZ equation in 1 + 1 and 2 + 1 dimensions with an Euler 
discretization scheme. The 1 + 1-dimensional case confirms the appearance of instabilities 
for high nonlinearities and large box sizes. These instabilities are suppressed by replacement 
of (Vh) 2 by a exponentially decreasing function of this quantity in the KPZ equation and 
subsequent discretization. Moreover, this change leads to consistent estimates of scaling 
amplitudes, in contrast to the discretization scheme of the original equation. In 2 + 1 
dimensions, we spanned a reasonable range of the model parameters where crossover effects 
(i. e. transients with EW growth or random growth) are not observed. We confirmed the 

10 



universality of roughness exponents, height distributions and roughness distributions in the 
steady state, which were previously obtained in discrete models. Estimates of skewness and 
kurtosis of those distributions were provided in order to show the quantitative agreement 
with previous results. We also showed evidence that the tails of the roughness distributions 
are stretched exponentials, which also agrees with discrete model results and suggests the 
non-Gaussianity of the steady state KPZ interfaces in 2 + 1 dimensions. 

Our results are not able to improve the accuracy of simulations of discrete models, which 
is expected for the computational limitations in the work with floating point operations. 
At first sight, this could seem to diminish the relevance of the present work. However, 
we stress that the connections between the lattice models and the KPZ equation do not 
follow from rigorous mathematical proofs, and some models are controversial at this point 
(e. g. ballistic deposition 



42j). Moreover, works on discrete models are not able to vary the 
parameters of the corresponding KPZ equations in a systematic way; indeed, only two or 
three of those models provide finite-size data which are clearly consistent with universality 



(see e. g. Ref. [13]). Thus, confirming results in the context of the KPZ equation itself, 
with different parameter values, is of great importance. As far as we know, this is the first 
quantitative discussion on height and roughness distributions obtained from integration of 
the KPZ equation in 2 + 1 dimensions. Those quantities are very useful for a complete 
characterization of a growth class, particularly due to the effects of scaling corrections in 
the estimates of exponents. 

We also believe that this work can motivate future studies in higher dimensions, where 
the debate on the existence of a finite upper critical dimension (the dimension where the 
nonlinear ity is always irrelevant) still remains Q, 13], despite the strong numerical evidence 



against it provided by results of two lattice models [16|, 1?|]. Since such studies would 
demand an efficient integration scheme to provide the best possible accuracy in reasonable 
simulation times, one must consider the possible advantages of other approaches, such as 



the pseudospectral methods [29|, [31] 
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TABLE I: Nonlinear parameters used in the integration of the KPZ equation in 2 + 1 dimensions 
and the corresponding constant for controlling instabilities. In all cases, v = 0.5 and a = 0.1. 

Set g c 
A 12 0.1 
B 24 0.5 
C 48 1.0 
D 96 4.0 
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FIG. 1: Interface profiles in d = 1 obtained in the integration of the KPZ equation with Eq. 
in times (a) t = 55.5 and (b) t = 55.8. Notice the different vertical scale in (a) and (b) due to the 
rapid growth of an instability in x ~ 120. 
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FIG. 2: Amplitude of squared roughness as a function of inverse box size obtained in the integration 



of the KPZ equation in d = 1 with Eq. 
/ ((V/i) 2 ) with c = 1 (squares). 



(circles) and modified Eq. ([8]) replacing (V/i) 2 by 
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FIG. 3: Effective roughness exponents as a function of inverse box size for KPZ interfaces in 
2 + 1 dimensions. The symbols correspond to sets A (diamonds), B (squares), C (circles), and 
D (crosses). For L < 64, error bars are smaller than the size of the data points. For L = 128, 
uncertainties in a(L) are near 0.01. 
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FIG. 4: (a) Skewness and (b) kurtosis of the height distributions of KPZ interfaces in 2 + 1 
dimensions versus inverse box size. Each symbol corresponds to the same set of Fig. 3. Error bars 
in S are of the order of the size of the data points. Uncertainties in Q are near 0.02 for L < 64 
and near 0.05 for L = 128. 




5 10 

(w 2 -<w 2 >)/a 

FIG. 5: Scaled roughness distribution of KPZ interfaces in 2 + 1 dimensions for sets C and D. Each 
symbol corresponds to the same set of Fig. 3. 
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FIG. 6: Effective exponents j(x) versus 1/x 2 of roughness distributions of KPZ interfaces in 2 + 1 
dimensions. Each symbol corresponds to the same set of Fig. 3. 
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